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Abstract 

Germanium detectors have very good capabilities for the investigation of rare phenomena like the neutrinoless double beta decay. 
Rejection of the background entangling the expected signal is one primary goal in this kind of experiments. Here, the attainable 
background reduction in the energy region where the neutrinoless double beta decay signal of 76 Ge is expected to appear has 
been evaluated for experiments using germanium detectors, taking into consideration different strategies like the granularity of the 
detector system, the segmentation of each individual germanium detector and the application of Pulse Shape Analysis techniques 
to discriminate signal from background events. Detection efficiency to the signal is affected by background rejection techniques, 
and therefore it has been estimated for each of the background rejection scenarios considered. Finally, conditions regarding crystal 
mass, radiopurity, exposure to cosmic rays, shielding and rejection capabilities are discussed with the aim to achieve a background 
level of ~10 -3 c keV -1 kg -1 y _1 in the region of interest, which would allow to explore neutrino effective masses around ~40 meV. 

Key words: double beta decay, germanium detectors, granularity, segmentation, pulse shape analysis, sensitivity 
PACS: 14.60.Pq, 23.40.-s, 24.10.Lx, 29.40.-n 



1. Introduction 

Investigation of neutrinoless Double Beta Decay (DBD) 
can shed light on interesting pending questions like the ab- 
J-h solute values of the neutrino mass and the properties of 
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neutrinos under CP-conjugation [1]. Germanium detectors 
enriched in double beta decay emitter 76 Ge, offer important 
advantages for this investigation in comparison with other 
kind of detectors and nuclei [2,3]: excellent energy resolu- 
tion, high purity materials and powerful background rejec- 
tion capabilities, well established detector technologies, fa- 
vorable nuclear matrix element, high transition energy Q 
around 2039 keV [4], .... Indeed, germanium double beta 
decay experiments (Heidelberg-Moscow and IGEX) have 
provided the most restrictive bounds on the effective neu- 
trinos mass [5,6]. 

Both Heidelberg-Moscow and IGEX experiments, 
presently finished, followed similar strategies using several 
massive high-purity germanium detectors (HPGe) (around 
2 kg each) with close-end coaxial geometries together 
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with active and passive shieldings in deep underground 
laboratories. As a continuation of these experiments, next- 
generation projects as Majorana [7] and GERDA [8], have 
been proposed incorporating different innovations like us- 
ing segmented germanium detectors or operating naked 
crystals in cryogenic liquids, like nitrogen or argon. 

Important advances in high-purity germanium detector 
technologies have been achieved in the last years [9] , allow- 
ing the construction of large efficiency HPGe crystals and 
developing the monolithic segmentation technique which 
provides both interaction position and energy information. 
The analysis of pulse shapes in highly segmented germa- 
nium detectors (taking into account not only net signals in 
a segment but also induced transient signals in the neigh- 
boring segments) has revealed as a promising technique for 
three-dimensional position determination [10]. First results 
of the operation of segmented germanium crystals within 
DBD projects have recently been presented [11,12]. Impor- 
tant achievements have been obtained in other contexts for 
the spatial resolution of events (see for instance Refs. [13]- 
[16]), greatly surpassing the position determination capa- 
bilities of previous DBD germanium experiments based on 
pulse shape discrimination in conventional germanium de- 
tectors [17]. The main difference between real double beta 
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decay and background events is that the former leave gener- 
ally only one energy deposit inside the crystal (" monosite" 
events), whereas the latter leave in many cases more than 
one ("multisite" events). For this reason, thanks to the 
analysis of pulse shapes, only background events produc- 
ing one energy deposit in the crystal can be mistaken with 
real double beta decay events. One goal of this work is 
therefore to evaluate the background reduction attainable 
at present for germanium detectors in the region between 
2 and 2.1 MeV where the neutrinoless double beta decay 
signal of 76 Ge is expected to appear. Reduction of back- 
ground is based on the granularity of the detector system, 
on the segmentation of each individual germanium detec- 
tor and on the application of Pulse Shape Analysis (PSA) 
techniques to discriminate signal from background events. 

In underground experiments searching for rare events 
like the nuclear double beta decay, the background entan- 
gling the expected signal comes mainly from environmen- 
tal gamma radiations (including radon emissions) and neu- 
trons at the laboratory, radioactive impurities (cither pri- 
mordial or cosmogenically induced) in the materials of the 
experimental set-up and cosmic muons arriving even deep 
underground [18,19]. To explore effective neutrino masses 
around 40 meV, a background level ^10~ 3 c keV -1 kg -1 
y _1 in the region of interest must be achieved in germanium 
double beta decay experiments. The two main sources of 
the background registered in the region of interest for the 
neutrinoless double beta decay of 76 Ge have been identi- 
fied in previous experiments like IGEX to be cosmogenic 
activation of germanium detectors (mainly 68 Ge and 60 Co) 
[20] and external gamma background above 2 MeV coming 
from 232 Th and 238 U chains. In fact, in present projects of 
new experiments these two sources continue to be a domi- 
nant background [7,8]. Therefore, in this work the particu- 
lar background components thought to be the most signifi- 
cant ones will be taken into account: cosmogenic 68 Ge and 
60 Co produced in the germanium crystal and the 2614.5 
keV emission from 208 T1 in the 232 Th chain. Precise con- 
ditions necessary to achieve a background level ~T0~ 3 c 
keV -1 kg^ 1 y _1 in the region of interest from these sources 
will be discussed. Recently, a Monte Carlo study of the 
background achievable in the GERDA experiment by anti- 
coincidence cuts between crystals and segments has been 
published considering the main radioactive impurities in 
the set-up [21]; it has also been shown in the context of 
this experiment that muon-induced contribution to back- 
ground can be of ^10 -4 c keV -1 kg -1 y _1 provided that 
muon veto system is used [22] . 

The paper is organized in the following way. In Sec. 2 
the raw backgrounds expected from the relevant sources 
are studied, including results from a recent evaluation of 
the production rates of the most relevant cosmogenic prod- 
ucts in germanium detectors at sea level. Then, the analy- 
sis of the effects of each one of the three reduction strate- 
gies considered (granularity, segmentation and PSA) for 
the different background sources will be presented in Sec. 
3. Detection efficiency to the neutrinoless DBD signal and 



corresponding sensitivity will be studied in Sec. 4 assum- 
ing the different background rejection scenarios. Finally, 
results and conclusions will be discussed in Sec. 5. 



2. Raw backgrounds 

An attempt has been made to estimate the raw contribu- 
tions to the detector counting rates coming from the main 
sources of background taken into account, before applying 
any technique for background reduction. 

For this estimate, and for the study of the different back- 
ground reduction strategies presented in the next section, 
a set of Monte Carlo simulations using GEANT4 [23] pack- 
age has been developed. As the first goal was focused on 
the crystal geometry, a natural germanium crystal without 
shielding has been defined as detector. Internal contami- 
nations of the crystal are emitted homogeneously inside it, 
whereas for external sources, the corresponding gammas 
were emitted homogeneously and isotropically from the sur- 
face of an external sphere. Standard GEANT4 models, in- 
cluding those specific for low energy, have been used for all 
the processes, isotopes decays, and particles simulated. For 
every simulation made, position and energy of each inter- 
action produced in an event have been registered. This in- 
formation has allowed us to make different analysis of the 
obtained data in the Region of Interest (Rol) between 2 
and 2.1 MeV. It has been simulated a number of events big 
enough to obtain a negligible statistical error (below 2% for 
all the studies made) . 



2.1. Cosmogenic radioactivity 

Long-lived radioactive nuclei induced by the exposure 
of the materials of the set-up to cosmic rays at sea level 
(during fabrication, transport and storage) may become 
very problematic for rare event experiments which oper- 
ate in deep underground locations, using active and pas- 
sive shields and selecting carefully radiopure components. 
Therefore, materials must be kept shielded against the 
hadronic component of the cosmic rays to prevent cosmo- 
genic activation, flying must be avoided and the exposure 
on the surface of the Earth should be reduced as much as 
possible. Since these requirements usually complicate the 
preparation of experiments (for example, crystal growth 
and mounting of detectors) it would be desirable to have 
reliable tools to quantify the real danger of exposing the 
materials to cosmic rays. 

For these reasons, cosmogenic activation in germanium 
double beta decay experiments was specifically studied in 
Ref. [24] . Using these and previous results from the litera- 
ture, together with specific simulations of the response to 
the cosmogenic background of germanium detectors, the 
expected counting rate due to this effect in the neutrinoless 
DBD region has been evaluated. 
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Fig. 1. Comparison of excitation functions for 68 Ge (top) and for 
60 Co (bottom) in natural germanium by nucleons from different 
sources: measurements (Horiguchi'83 [25], Aleksandrov'91 [26] and 
Norman'05 [27]), calculations using the YIELDX code [28], and 
Monte Carlos results (elaborated by the Majorana Collaboration us- 
ing ISABEL [7] , taken from the MENDL2 libraries based on ALICE 
code [29] and taken from the library in Refs. [30,31] using CEM95 
and HMS-ALICE codes). 

2.1.1. Production rates 

Excitation functions of 60 Co and 68 Ge, as well as of 
other long-lived products induced in germanium, have been 
compiled using available measurements and different cal- 
culations taken from libraries or made on purpose in Ref. 
[24] (see Figure 1). In an attempt to find the most re- 
liable selection of the excitation functions, those calcula- 
tions with minimum deviations with respect to experimen- 
tal data were taken into consideration: below 150 MeV, 
HMS-ALICE results for neutrons [30,31], and above this 
energy, YIELDX calculations. 

Once chosen the excitation functions, production rates 
of 60 Co and 68 Ge were calculated using the cosmic neutron 
spectrum from Ziegler [32], for both natural and enriched 
germanium (86% of 76 Ge, 14% of 74 Ge). Results are sum- 
marized in Tables 1 and 2, compared with some previous 
estimates. 

A production rate of ~5 kg~ 1 d~ 1 can be safely consid- 
ered in natural Ge for 60 Co while for 68 Ge the production 
rate of ~90 kg _1 d _1 found in this work is significantly 
higher than in previous estimates including the measure- 
ments in Ref. [35] due to the large contribution of neutrons 
below 100 MeV. Since there is now no available measure- 
ment of production cross sections by neutrons for the rele- 
vant isotopes, it is difficult to assess to what extent the two 
basic assumptions of the presented calculations are valid. 



For enriched detectors, production of 68 Ge, although much 
more reduced than the one for natural germanium, has been 
found one order of magnitude higher than previous esti- 
mates commonly used [34,35]. 60 Co production seems to be 
similar in enriched and natural material. 

2.1.2. Counting rates in the region of interest 

A set of Monte Carlo simulations based on GEANT4 code 
have been made to reproduce the response of cylindrical 
germanium detectors of different masses to the background 
sources, including 60 Co and 68 Ge isotopes induced in the 
germanium crystal. 

The number of events registered in the 2-2.1 MeV Rol per 
isotope decay can be deduced using these MC simulations. 
For 2(4)-kg detectors, these numbers are 0.00016 (0.00020) 
c/keV/decay of 60 Co and 0.00023 (0.00029) c/keV/decay 
of 68 Ge (in agreement with the ones deduced in simulations 
made by the Majorana [7] and GERDA [8] collaborations). 

Using this information from GEANT4 simulations to- 
gether with the production rates previously presented, the 
counting rates in the region of interest due to the cosmo- 
genic contaminations can be derived for certain exposure 
and cooling times of the material. To properly compare re- 
sults, the same times typically used in the GERDA project 
[36] will be considered: 30 days of exposure to cosmics rays 
for 60 Co production and 180 days of exposure plus 180 days 
of cooling for 68 Ge. Table 3 presents the obtained counting 
rates in these conditions and using production rates of 5 
kg _1 d _1 for 60 Co and 1 kg _1 d _1 for 68 Ge; it is also shown 
the estimate when considering a 68 Ge production rate of 10 
kg _1 d _1 and 2 years of cooling. As it can be seen, for this 
high 68 Ge production rate it would be necessary to wait 
more than two years (instead of 3 months) to achieve count- 
ing rates of the same order than when considering the low 
production rate. Rates presented in Table 3 correspond to 
conventional germanium detectors, that is, without taking 
into account neither segmentation nor PSA techniques. 

2.2. 2614.5 keV 20S Tl emissions 

An evaluation of the expected background coming from 
the 2614.5 keV line from 208 T1 has been made, consid- 
ering both the environmental gamma background in the 
laboratory and the 232 Th intrinsic radioimpurities in the 
lead shielding expected to be surrounding detectors. The 
response to 2614.5 keV photons for the counting rates in 
the 2-2.1 MeV Rol in 2-kg and 4-kg germanium detectors 
has been estimated by Monte Carlo simulation. A flux of 
~0.1 cm~ 2 s _1 for environmental 2614.5 keV photons has 
been assumed, according to recent measurements in the 
new Canfranc Underground Laboratory [37], and an activ- 
ity of 1 /iBq/kg of 232 Th in lead has been considered just as 
a reference value. In these estimates, a spherical cavity with 
radius R=30 cm for placing detectors inside has been sup- 
posed, and two different shielding configurations with 30 



Table 1 

Production rates (in kg d -1 ) in natural germanium obtained considering HMS-ALICE and YIELDX below and above 150 MeV respectively, 
together with previous estimates. Some results from Ref. [35] are based on Monte Carlo calculations (MC) while others come from measurements 
(exp). 





HMS-ALICE+ YIELDX [24] Ref. [33] Ref. [34] Ref. [35] (MC) Ref. [35] (exp) 


68 Ge 
60 Co 


77+12=89 58.4 26.5 29.6 30±7 
0.3+4.5=4.8 6.6 4.8 



Table 2 

As Table 1, but for enriched germanium (86% of 76 Gc and 14% of 74 Gc). 





HMS-ALICE+ YIELDX [24] Ref. [34] Ref. [35] 


68 Ge 
60 Co 


2.8+10=13 1.2 0.94 
0.02+6.7=6.7 3.5 



Table 3 

Estimates of counting rates R (in units of 10~ 3 c kcV -1 kg _1 y _1 ) in the 2-2.1 MeV Rol from cosmogenic contaminations in 2-kg and 4-kg 
germanium detectors. Uncertainties in these estimates arc discussed in the text. 
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Table 4 

Estimates of counting rates (in units of 10~ 3 c keV -1 kg _1 y _1 ) in the 2-2.1 MeV Rol from 2614.5 keV 208 T1 emissions in 2-kg and 4-kg 
germanium detectors (see text). Two shielding configurations with 30 and 40 cm of lead have been taken into consideration. Uncertainties 
in these estimates arc discussed in the text. 
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2.2 



and 40 cm of lead surrounding the cavity have been taken 
into account. Table 4 summarizes the obtained results. 

For the intrinsic radioimpurities in the lead shielding, 
it has been checked that around 90% of the events come 
from the most internal 5 cm of shield, and in fact, virtually 
the same counting rates have been found when considering 
30 or 40 cm of lead (despite the significant difference in 
mass and consequently in 232 Th activity between the two 
configurations). For the external gamma background, the 
additional 10 cm when assuming a 40-cm-thick lead layer 
reduce the counting rates around two orders of magnitude. 
Comparing the background levels expected in 2 and 4-kg 
detectors in Tables 3 and 4, the largest crystals have a lower 
background level (~20%) from 208 T1 but higher counting 
rates (~30%) from internal cosmogenic impurities. This 
fact will be discussed at Sec. 3.1. 

Uncertainties in the estimates of raw backgrounds pre- 
sented in Tables 3 and 4 come in principle from the sim- 
ulations of the response of detectors to the different back- 
ground sources as well as from the inputs considered for 
production rates of cosmogenic isotopes and levels of 208 T1 
photons. Counting rates are directly proportional to these 
inputs, which are thought to be the main source of error 
since GEANT4 can reproduce electromagnetic processes 



with a few per cent error and statistical errors in simulations 
have always been kept below 2% as stated before. Produc- 
tion rate of 60 Co has an uncertainty of ~50% (when con- 
sidering different available estimates) and for 68 Ge it could 
be up to one order of magnitude. Environmental gamma 
fluxes depend on the particular underground location. For 
232 Th impurities in very pure lead just upper bounds have 
been derived (see Refs. [38,39]) and therefore the assumed 
value must be taken just as a reference. 



3. Rejection of background events 

The main goal in the study of the different ways to re- 
duce the background was to try to analyze the correlation 
between the experiment design and background level in the 
region of interest. Assuming a hypothetical final set up of 
an experiment with a total mass of germanium around some 
tens of kilograms, the aim was to determine what is the op- 
timal mass distribution and features of detectors to have a 
background level as low as possible. Three topics that can 
determine the best configuration, as we pointed out before, 
are: granularity of the experiment detectors, segmentation 
of the crystals and analysis of the obtained pulses. 
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3.1. Granularity 

The first step to determine the best set up of the experi- 
ment consists in analyzing what could be the optimal mass 
of the detectors that build the whole experiment to have 
a background level as low as possible, if a fixed total mass 
of germanium is assumed. With this purpose, cylindrical 
crystals with the same value for diameter and height and 
masses between 0.1 and 4 kg have been simulated. 

If only the events with an energy deposit between 2 
and 2.1 MeV are taken into account, we can see how for 
internal contaminations and a given specific activity, the 
higher mass detectors register a higher background level 
(see Figs. 2 a,b). 4- kg detectors register 28% (26%) more 
events for 60 Co ( 68 Ge) than 2-kg ones. This is not the case 
for external contaminations, because for a given activity, 
the heavier detectors have a lower background (see Fig. 2 
c). 4-kg detectors register 54% less events coming from 
2614.5 keV photons than 2-kg detectors. These values con- 
firm that the optimal configuration of an experiment de- 
pends on what kind of background we want to reduce more, 
that coming from internal contaminations or from external 
ones, taking also into account that the dependency between 
the mass of the detectors and the background events reg- 
istered is stronger for the external contamination. These 
data together with conclusions obtained from the study of 
the segmentation of the crystal and pulse analysis, as ex- 
plained later, will determine the best configuration. 

3.2. Segmentation 

Regardless of the mass of the used detectors, other way to 
reduce the background level is by segmentation of the crys- 
tals and further application of anticoincidence techniques 
between segments. The aim was to quantify the maximum 
background reduction that could be obtained from segmen- 
tation of the crystals. For the cylindrical detectors studied, 
two different ways to divide it were considered: segmen- 
tation in planes and segmentation in sectors (longitudinal 
and transversal segmentation respectively) (see Fig. 3). 

Table 5 shows how the reduction of counts registered be- 
tween 2 and 2.1 MeV in detectors of 2 and 4-kg is bigger 
for higher number of segments, combining longitudinal and 
transversal segmentation. For a 4-kg detector with the high- 
est number of segments studied (66 segments distributed in 
11 transversal slices and 6 angular sectors), approximately 
only 2 out of 100 events will not be rejected by anticoin- 
cidence techniques for 60 Co contamination. For 68 Ge the 
ratio is around 5 out of 100 events and for 2614.5 keV ex- 
ternal gammas, less than a half of total events. This reduc- 
tion can be observed qualitatively in the spectra obtained 
from simulations for the different background sources stud- 
ied (see Fig. 4). 

It is clear that for higher number of segments in the crys- 
tal, we are able to reject more background events. But the 
segmentation of a crystal is limited for some reasons. First 
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Fig. 2. Background level in the 2-2.1 MeV Rol depending on the 
component detector mass. For internal contaminations, 60 Co (a) and 
68 Ge (b), it is represented in counts per keV every 10 4 decays per kg. 
For external photons of 2614.5 keV (c), it is represented in counts 
per keV per kg every 10 4 interacting photons. 



of them is the difficulty associated to the reduction of the 
width of the transversal segments. It is also necessary to 
make a segmentation to reject a high number of background 
events without losing too much efficiency in the double beta 
decay events detection. The different segmentation config- 
urations studied have transversal segments with a width 
of 1 cm approximately. This size is reachable using actual 
segmentation techniques and provides a good efficiency for 
the detection of double beta decay events, as it will be dis- 
cussed later in Sec. 4.1. 

Considering the biggest segmentation of all studied and 
applying the reduction factors obtained and shown in Ta- 
ble 5 to the raw rates for different background sources esti- 



Table 5 

Percentage of rejected events by crystal segmentation and anticoincidence techniques for registered events between 2 and 2.1 MeV. Different 
segmentation schemes for 2 and 4-kg detectors have been considered. 
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Fig. 3. Segmentation scheme for a germanium crystal with 9 transver- 
sal slices and 6 longitudinal sectors to obtain 54 segments. 

mated in Sec. 2, the resulting rates can be derived and are 
presented in the first and third columns of Table 6. 



3.3. PSA 

Besides the improvements made in the detector, like the 
segmentation of the crystal, a reduction of the background 
level from the analysis of the pulses registered can also 
be obtained. For this purpose, as pointed out before, it 
is necessary to distinguish between background and real 
double beta decay events in order to reject the first ones. 
Application of PSA in segmented crystals allows to increase 
the spatial resolution of conventional germanium detectors 
to obtain the position of energy deposits with a very good 
accuracy in all the dimensions [16] or at least a correct 
identification of the number of interaction points [40] . 

The data obtained in the simulations can be reanalyzed 
for a given spatial resolution, grouping all the partial energy 
deposits with a separation lower than the resolution and 
considering these groups like indivisible energy deposits. 
Then, it is possible to determine how many of these indivis- 
ible deposits each background event has. In figure 5, distri- 
butions of the number of energy deposits per event for dif- 
ferent background sources and 4-kg detectors are showed, 
assuming 2 different values for the spatial resolution, 3 and 
5 mm, which seem to be at reach today according to the 
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Fig. 4. Comparison between the spectra obtained for the different 
contaminations studied: 6t) Co (a), 68 Ge (b) and 2614.5 keV photons 
(c), for a 4-kg detector without segmentation (black line), with 11 
transversal segments (blue line) and with 66 segments, 11 transversal 
by 6 longitudinal (red line). 



Table 6 

Estimates of counting rates (c keV - 1 kg - 1 y — 1 ) in the 2-2.1 MeV Rol from cosmogenic contaminations and 2614.5 keV gamma emissions 
from 208 Tl in 2-kg and 4-kg germanium detectors, assuming maximum reduction factors deduced for segmented crystals using anticoincidence 
techniques and for segmented crystals using PSA. Same conditions for cosmogenic production and 2614.5 keV emissions than in Tables 3 
and 4 are taken into account. 
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work developed in Ref. [160; it is worth noting that these 
distributions have been obtained without making any par- 
ticular definition of the segmentation scheme. The worst 
value of the spatial resolution can give an idea of the loss 
of rejection efficiency depending on the spatial resolution 
finally achieved. In Table 7, rejection factors after elim- 
inating "multisite" events, between 2 and 2.1 MeV, are 
showed in function of the 2 values considered for the spa- 
tial resolution and for all background sources studied and 
2 and 4-kg detectors. For the heaviest crystal and assum- 
ing a 3 mm spatial resolution, it is possible to reject 99.9% 
of background events coming from 60 Co, 99.2% from 68 Ge 
and 60.3% from 2614.5 keV photons; the values are 99.5%, 
97.8% and 56.5% respectively if a resolution of 5 mm is 
considered. 

Spatial resolution of the detectors allows to reject more 
or less background events depending on how the energy of 
these events is distributed in the crystal. One way to predict 
what could be the maximum reduction factor consists in 
determining the maximum distance between all the energy 
deposits of the same event (that we call maximum interdis- 
tance D max , see Eq. 1) and studying how it is distributed. 
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(1) 



This maximum interdistance depends on different factors 
like the mass of the detector, the scheme of the decay of the 
isotope that produces the background event or the origin 
of the event, because it is different if the contamination is 
located inside the crystal or if comes from outside of the ex- 
perimental setup. In Fig. 6, the distribution of these max- 
imum interdistances for 2 and 4-kg detectors and all the 
background sources studied are presented, showing how all 
the factors mentioned previously have influence on these 
distributions. From these plots in Fig. 6, rejection factors 



1 A genetic algorithm for the decomposition of multiple hit events is 
presented in Ref. [16], considering the features of a cylindrical closed- 
end germanium detector with a mass of almost 2 kg and 6 angular 
sectors and 4 transversal slices. However, it is reported that the 
approach has no limitation concerning the geometry of the crystal, 
the number and layout of the segments or the number of interactions. 



for any considered experimental spatial resolution can be 
deduced knowing that all events with a maximum interdis- 
tance lower than the spatial resolution will be labelled as 
"monosite" like double beta events. 

Second and fourth columns in Table 6 show the expected 
counting rates in the 2-2.1 MeV Rol for background con- 
tributions in 2-kg and 4-kg germanium detectors, applying 
the reduction factors obtained when considering a 3D spa- 
tial resolution of 3 mm from PSA on the raw backgrounds 
estimates in Sec. 2. Total background levels have been cal- 
culated (see last row in Table 6) including cosmogenic 60 Co 
and 68 Ge in the detector as well as external and lead 2614.5 
keV emissions and using the most favorable conditions, that 
is, adding values in first, second, fifth and sixth rows. 

It must be noted that no technical limitation in the 
process of rejecting events thanks to coincidences between 
segments or PSA with a certain spatial resolution has been 
taken into account. Consequently, the rejection factors 
summarized in Tables 5 and 7 must be considered as the 
best ones achievable and give the limit of the power of these 
rejection techniques. Uncertainties in the total background 
levels presented in Table 6 must be dominated by the un- 
certainties in the raw backgrounds, discussed in Sec. 2, 
since statistical errors in the simulations are much smaller. 

4. Efficiency to neutrinoless DBD signal and 
sensitivity 

The sensitivity of a neutrinoless DBD experiment is often 
evaluated using the detector factor-of-merit F^> defined as: 



F D = 4.17 x 10 



26 



/ MT 



W V bAE 



(2) 



with f the isotopic abundance of the DBD emitter, W the 
atomic weight of the source material, MT the exposure of 
the experiment, b the background level (typically expressed 
in counts per keV, per kg and per y) , AE the energy win- 
dow where the signal is expected to appear, dependant on 
the energy resolution of the detector, and e the detection 
efficiency. Fn is interpreted as the lifetime corresponding 
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Table 7 

Percentage of rejected events in the 2-2.1 MeV Rol by PSA techniques considering a spatial resolution of 3 and 5 mm for 2 and 4-kg detectors. 



2kg 

3 mm 5 mm 
resolution resolution 



60 Co 


99.7 


99.0 


99.9 


99.5 


68 Ge 


98.9 


97.0 


99.2 


97.8 


external 2614.5 keV gammas 


55.6 


51.2 


60.3 


56.5 



4kg 

3 mm 5 mm 
resolution resolution 
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Fig. 5. Distribution of the number of energy deposits per event in the 
2-2.1 MeV Rol for all the background contributions studied: 60 Co 
(a), 68 Ge (b) and 2614.5 keV photons (c), for 4-kg detector and 
considering a spatial resolution of 3 (solid line) and 5 mm (dashed 
line). Monosite events are singled out in red (blue) for 3 (5) mm 
resolution. 
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Fig. 6. Distribution of the maximum interdistance between energy 
deposits of the same event in the 2-2.1 MeV Rol for 60 Co (a)and 68 Ge 
(b) internal contaminations and for 2614.5 keV external photons (c), 
for 2-kg (red line) and 4-kg (black line) detectors. 

to the minimum detectable number of events over a back- 
ground at 1(7 confidence level. 

The neutrino effective masses which can be explored by 
an experiment with a detector factor-of-merit Fd can be 



determined as: 



< m v > z 



FnF; 



D r N 



(3) 



with m e the electron mass and Fjv the nuclear factor-of- 
merit, defined as Fjv=G 0,y |M 0,y | 2 , being G°" a kinematical 
factor and M 0l/ the nuclear matrix element qualifying the 
likeliness of the transition. 

The background level b achievable in germanium experi- 
ments has been analyzed in previous sections under differ- 
ent background reduction schemes; but the application of 
anticoincidence rejection or PSA techniques unfortunately 
affects also the efficiency for the detection of the neutrino- 
less DBD signal. Therefore, a study of this efficiency has 
been carried out and is presented here. 

4.1. Efficiency to signal 

In order to evaluate the dependency between background 
reduction techniques and loss of efficiency, it is necessary to 
simulate neutrinoless DBD (0vf3f3) events inside the detec- 
tor to apply them the same treatment that a background 
event. This procedure allows to estimate the percentage 
of 0vf3(3 events rejected losing detection efficiency. Signal 
events can be missed either because of the escape of the 
Brcmsstrahlung radiation produced by electrons or because 
0vf3f3 events are mistaken as background due to the spatial 
distribution of their energy deposits. 

To define Qv(3(3 events in the simulations, two electrons 
with sum energy of 2040 keV (Q value for 76 Ge DBD) have 
been considered, neglecting in first approach the angular 
correlation between electrons. In order to check if this an- 
gular correlation could have substantial influence to detect 
a 0v(3(3 event, a first study was made simulating two elec- 
trons of 1020 keV each emitted randomly, with the same 
and with opposite directions. Decays are distributed ho- 
mogenously in the crystals, for 2 and 4-kg detectors. Tables 
8 and 9 show that the differences between the detection ef- 
ficiency factors in the cases described before are less than 
1% for the different angular correlations considered, vali- 
dating the approximation taken. 

Another point in the simulation of 0is(3/3 events is the en- 
ergy distribution of the two electrons. Four different config- 
urations were studied: two electrons of 1020 keV each (half 
of the total energy), two electrons of 1500 and 540 keV, two 
electrons of 1734 and 306 keV (75% and 25% of the full en- 
ergy respectively) and one electron of 2040 keV. It is impor- 
tant to point that the last case is not a real case but could 
be useful to obtain limit values in the detection efficiency 
because Bremsstrahlung probability is higher for more en- 
ergetic electrons. For all these energy configurations, the 
electrons were emitted in random directions using the same 
simulation package that for the background events. Figure 
7 shows an example of the spectrum registered in the de- 
tector after the simulation of Of/3/3 events. In this case in 
particular, two electrons of 1020 keV each were emitted in 
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Fig. 7. Energy spectrum registered for simulated 0^/3/3 events com- 
posed by 2 electrons of 1020 keV each emitted in random directions 
in a 4-kg detector. 



a 4-kg detector. Three different regions can be identified in 
this spectrum: 

- Up to 1020 keV, energy is lost from both electrons, either 
by escape of Bremsstrahlung radiation and/or escape of 
the electron itself (if the decay is close enough to the 
detector surface). The probability that this happens for 
both electrons is quite low and for this reason a step 
appears at 1020 keV. 

- From 1020 keV up to the peak, events where one elec- 
tron deposits all the energy and the other one suffers 
Bremsstrahlung losses or escapes from the detector, are 
registered. In a smaller percentage, events where both 
electrons do not deposit the full energy can appear here 
too. 

- In the peak, electrons have deposited the full energy. 
This is the region where the 0vf3f3 signal is expected. As 
explained later, some events in the peak can be labelled as 
"multisite" events and therefore rejected as background, 
when electrons produce Bremsstrahlung emission which 
is finally absorbed in the crystal. 

Using the simulation of Qv(3[3 events, the relationship be- 
tween the detection efficiency and the detector mass, di- 
rectly related with the granularity of a future experiment, 
can be analyzed. Tables 8 and 9 show, for two 1020 keV elec- 
trons, that the difference between the detection efficiency 
for 2 and 4-kg detectors is less than 1% for all the back- 
ground rejection configurations, but always better for the 
heaviest detector. For this reason, the studies of the differ- 
ent energy configurations were made for a 4-kg crystal. 

Table 10 summarizes the detection efficiencies obtained 
in a 4-kg detector for the different energy configurations in 
Qv(3f3 events. Two main conclusions can be obtained from 
this table. One is that PSA techniques offer better detection 
efficiency than background rejection based just on crystal 
segmentation. The explanation is that the anticoincidence 
between segments rejects all the 0vf3f3 events with energy 
deposits in two or more segments of the crystal, while PSA 
rejects events with two separated energy deposits, typi- 
cally due to Brcmsstrahlung but, in principle, could allow 
events close to the borders. A useful information to under- 
stand this is the maximum distance between all the energy 
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Fig. 8. Distribution of the maximum interdistance between energy 
deposits of the same 0^/3/3 event at the peak in a 4-kg detector for 2 
electrons of 1020 keV each (green), 1500 + 540 keV electrons (red), 
1734 + 306 keV electrons (blue) and one 2040 keV electron (black). 



Fig. 10. Exposure MT (in kg-y) necessary to reach the sensitivity 
to explore neutrino effective masses < m„ > of 40 meV, versus the 
detection efficiency to Ou/3/3 events considering full crystals (solid 
line), anticoincidences in 11 x 6 segmented crystals (dashed line), and 
PSA with a spatial resolution for PSA of 3 mm (dot and dash line). 
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Fig. 9. Detection efficiency versus the energy of the most energetic 
electron of a Ov/3/3 event in a 4-kg detector considering a full crystal 
detector (black points), anticoincidences in 11 x 6 segmented crystals 
(blue squares), and PSA with 3 mm (red triangles) and 5 mm (white 
circles) spatial resolution. 

deposits obtained in the simulation for Qu(3f3 events (the 
maximum interdistance defined in Section 3.3 and calcu- 
lated using Eq. 1). Figure 8 shows the distribution of the 
interdistances for signal events and for all the energy con- 
figurations studied in a 4-kg detector. In all the cases, most 
of the events have a maximum interdistance below 3 mm, 
ensuring that they will be considered as " monosite" events 
by PSA and can be separated safely from background ones. 

The other important point from Table 10 is the depen- 
dency between the efficiency and the energy of the most 
energetic electron of the Qv(3j3 event. The more equal the 
electron energies, the higher the detection efficiency. This 
relationship is logical due to the Bremsstrahlung probabil- 
ity, that grows proportional to the energy of the electron. 
For this reason, the non real case of a 2040 keV electron 
can be useful to estimate the lower limit of the detection 
efficiency. The obtained dependencies of the detection ef- 
ficiency on electron energy can be fitted to a grade two 
polynomial, for the different background rejection config- 
urations considered (as shown in Fig. 9); convoluting this 
polynomials with the single electron spectrum of Quf3[3 de- 
cays in 76 Ge, an overall value for the efficiency to signal can 
be estimated. Last row in Table 10 presents these results. 



4.2. Sensitivity 

The overall detection efficiency to signal estimated above 
has been used together with Eqs. 2 and 3 to compare the 
sensitivity of germanium DBD experiments using different 
background rejection schemes. Three different situations 
have been taken into consideration: no rejection technique, 
the application of anticoincidence rejection for crystals with 
6x11 segments, and the use of PSA techniques assuming a 
3 mm spatial resolution. Detectors enriched in 76 Ge at 86% 
and total exposures of MT=100 and 1000 kg-y have been 
assumed. An energy window of 3.5 keV has been considered, 
as in Ref . [3] . Values of background level and corresponding 
detection efficiency deduced in this work for 4-kg crystals 
have been used for each situation. Table 11 summarizes 
the evaluated sensitivities, presenting the Yd values and 
the corresponding effective neutrino masses considering the 
average nuclear factor-of-merit Fjv=7.3x 10~ 14 y _1 used in 
Ref. [3]. Neutrino masses below 50 meV can be explored 
with very segmented crystals or applying PSA techniques 
and for high enough exposures. 

To achieve a certain sensitivity to the effective neutrino 
mass, the required exposure of an experiment depends on 
the background level in the region of interest and the de- 
tection efficiency to the signal (see Eqs. 2 and 3); both 
background level and efficiency are different when differ- 
ent background rejection schemes are considered. Figure 10 
shows this dependency between the exposure and the effi- 
ciency in the three background rejection schemes consid- 
ered before. Using the overall signal efficiency values pre- 
viously estimated, it can be deduced from this plot that 
despite the loss of efficiency, segmented crystals working 
in anticoincidence require an exposure one order of magni- 
tude lower than that of non segmented detectors to explore 
the same range of neutrino masses. If PSA techniques are 
used with 3 mm of energy resolution, an additional factor 
of ~two of reduction is achieved in the necessary exposure. 



Table 8 

Detection efficiency (%) for 0^/3/3 events composed by two electrons of 1020 keV each emitted in random, opposite and same directions. 
Different background rejection configurations arc considered, for 2-kg detectors. 
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Table 9 

As table 8, but for 4-kg detectors. 
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Tabic 10 

Detection efficiency (%) for Oi^/3/3 events composed by two electrons emitted in random directions with different energy schemes. Different 
background rejection configurations arc considered for 4-kg detectors. Last row shows an overall efficiency taking into account the single 
electron 0^/3/3 spectrum. 
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Tabic 11 

Comparison of the sensitivity of germanium DBD experiments under different background rejection schemes, evaluated following Eqs. (2) 
and (3) (see text). 
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5. Discussion and conclusions 

Some conclusions can be drawn from the analysis of the 
different strategies for background reduction in germanium 
double beta decay experiments. The study of the granular- 
ity of the detector system shows that heavier crystals are 
better to reduce the contribution of external radioimpuri- 
ties, but worse to reduce background coming from internal 
contaminations. By applying the most powerful segmen- 
tation techniques taken here into consideration in a 4-kg 
detector, 2(5) out of 100 events due to internal impurities 
from 60 Co ( 68 Ge) would remain in the Rol, while for ex- 
ternal contaminations, about half of the events would be 
rejected. A 3-dimensional spatial resolution of 3 mm, ob- 
tained by means of PSA in segmented detectors, would 
allow to reject more than 99% of background events due 
to cosmogenic isotopes induced in the crystal, and around 
60% of those coming from external 2614.5 keV photons. 

According to numbers presented in Table 6, a background 
level of 1.1 (1.6) 10" 3 c keV" 1 kg^y" 1 due to the stud- 
ied background sources could be achieved using 4 (2)-kg 
crystals when considering very precise conditions. For the 
production of 60 Co ( 68 Ge), the exposure time is of 30 (180) 
days and the production rate is 5 (1) kg _1 d _1 . In the case 
of 68 Ge, a cooling time of 180 days has been also taken 
into account. For the external 2614.5 keV, an environmen- 
tal flux of 0.1 cm~ 2 s _1 has been assumed and the use of 
a 40-cm-thick lead shield considered. For intrinsic 232 Th 
impurities in the lead shielding, the activity supposed is 1 
/iBq/kg. In these optimal conditions, the raw background 
is reduced by more than one order of magnitude thanks to a 
3 mm spatial resolution achieved by PSA in segmented de- 
tectors. The most relevant contribution is that of the 2614.5 
keV produced in the lead shielding itself by 232 Th impuri- 
ties, followed by the one from cosmogenic 68 Ge induced in 
the crystal. External 2614.5 keV photons and cosmogenic 
60 Co are almost negligible. The use of large crystals, hav- 
ing less background for external contaminations than the 
small ones, seems more adequate since background rejec- 
tion by anticoincidence between segments or by PSA is effi- 
cient enough in the reduction of contributions from internal 
radioimpurities . 

Some comments are in order if considering other less fa- 
vorable conditions for the raw background or the rejection 
strategies: 

- Contribution from 60 Co emissions is in general negligible 
thanks to the very good rejection factors attainable either 
with PSA in segmented detectors or just anticoincidence 
techniques between segments of the detectors. 

- 68 Ge contribution in the scenario described above con- 
sidering a production rate of 1 kg _1 y _1 is very small. 
However, it must be noticed that the production rate cal- 
culated in different estimates presented in Sec. 2.1.1 is 
much higher. It has been shown that if the production 
rate was 10 times bigger, a cooling time of 2 years instead 
of 3 months would be necessary to achieve background 



contributions of the same order of magnitude. 

- To shield external 2614.5 keV photons from environmen- 
tal gamma radiation in the laboratory, 40 cm of lead are 
mandatory. When using just 30 cm, the goal of a back- 
ground level of 10~ 3 c keV -1 kg _1 y _1 cannot be fulfilled. 

- Regarding intrinsic 232 Th impurities in the lead shield, 
it is worth noting that the achievement of a radiopurity 
of 1 /iBq/kg cannot be taken for granted since common 
upper limits to contaminations of this chain achieved in 
lead are of some hundreds of /iBq/kg (see for instance the 
ILIAS database on radiopurity measurement^ and even 
the best limits obtained with the most modern germa- 
nium spectrometers are of some tens of /iBq/kg [38,39]. 
Therefore, to further reduce the background level from 
this source, improvement of the radiopurity of this ma- 
terial or using radiopure copper instead of lead in the 
inner part of the shielding would be very useful. It has 
also been proved that increasing the thickness of the lead 
layer is not relevant for these shielding radioimpurities. 

- For the spatial resolution, if a worse performance of 5 
mm was achieved, the best total background level would 
be of 1.5 (2.0) 10" 3 c keV" 1 kg^d" 1 for 4 (2)-kg crys- 
tals. No dramatic difference in the background level from 
2614.5 keV would be produced, while an increase in a fac- 
tor of ~3 (4) would be registered for 68 Ge ( 60 Co) emis- 
sions, which, fortunately, seem not to be the dominant 
background source. 

- The use of 4-kg crystals gives a more reduced back- 
ground, but 2-kg detectors, in use for a long time, could 
be acceptable. Segmentation schemes assumed are feasi- 
ble for present germanium detector technologies and give 
a best total background level of 1.9 (2.2) 10~ 3 c keV" 1 
kg _1 y" 1 for 4 (2)-kg crystals. Therefore, the option of 
applying just segmentation techniques must not be dis- 
carded. 

The overall detection efficiency to neutrinoless DBD sig- 
nals has been evaluated by Monte Carlo simulation for the 
different background rejection scenarios, finding for 4-kg 
crystals a reduction from ^94% to ~80% when consider- 
ing anticoincidences in 11x6 segments and to ~85% if PSA 
techniques are applied with a 3 mm spatial resolution. 

The sensitivity of DBD experiments depends on both 
achieved background level in the Rol and detection effi- 
ciency; it has been shown that the improvement in the for- 
mer thanks to rejection techniques largely compensate the 
loss in the latter since experiments with these techniques 
require much lower exposure for a fixed sensitivity. 

In summary, it seems that contribution from dominant 
background sources in previous germanium double beta de- 
cay experiments could be reduced down to 10~ 3 c keV -1 
kg~ 1 y _1 in the Rol using present detector technologies, 
which allows to explore effective neutrino masses even be- 
low 50 meV. 



2 Available at http://radiopurity.in2p3.fr/ and developed within 
the ILIAS (Integrated Large Infrastructures for Astroparticle Sci- 
ence) EU project. 
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